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Abstract 



By using bosonization method and unitary transformation, we give a gen- 
eral relation between the dimensionless tunneling conductance and the frac- 
tional Coulomb blockade conductance peak splitting which is valid both for 
weak and strong transmission between two quantum dots, and show that the 
tunneling conductance has a linear temperature dependence in the low energy 
and low temperature limit. 
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Electron tunneling in a mesoscopic structure may be drastically influenced by the charg- 
ing effects. If the charge spreading is impeded by weak links, or by a special geometry of 
the structure, the charging suppresses the electron tunneling. Such a suppression of the 
electron tunneling is commonly called as the Coulomb blockade |1|;@,[J. It has become pos- 
sible to observe the Coulomb blockade effect in semiconductor heterostructures where the 
geometry of the system can be easily modified by adjusting the voltages on special gate 
electrodes. Recently, the electron tunneling through two quantum dots connected by single 
quantum point contact which is controlled by adjusting voltage on additional gate electrode 
has been extensively studied, both experimental @J^,|,0,H and theoretical P, [lO| . [TT| . [T2| , p!3| . 
The electron tunneling between two quantum dots leads to a decay of the Coulomb blockade 
of the individual dot. For a pair of electrostatically identical quantum dots, the progress of 
this decay can be described by tracking the splitting of the Coulomb blockade conductance 
peaks as they evolve from doubly degenerate single-dot conductance resonances to that of 
nondegenerate double-dot peaks with twice the original periodicity. The Coulomb block- 
ade peak splitting significantly depends upon the transmission coefficient T of the quantum 
point contact, i.e., the tunneling conductance. For the case of T = 1, the Coulomb blockade 
between two quantum dots disappears, they become a double-quantum-dot. For the case 
of T = 0, the Coulomb blockade is maximum which completely suppresses the electron 
tunneling between two quantum dots, therefore, they are separated from each other. 

In this paper, by using bosonization method and unitary transformation developed in 



Ref. ||14|| , we can effectively treat the system of two electrostaticaly identical quantum dots 
connected by one quantum point contact which is controlled by adjusting voltage on ad- 
ditional gate electrode, and give a general relation between the dimensionless tunneling 
conductance g and the fractional Coulomb blockade peak splitting /, which is valid both for 
weak (T ~ 0) and strong (T ~ 1) transmission coefficients. We also study the temperature 
dependence of the tunneling conductance in the low energy and low temperature limit. 

In generally, each individual quantum dot can be described by two-dimensional electron 
gas, and the Coulomb charging energy E c = e 2 /(2C), where C is the capacitance of the 



individual quantum dot, is very large compared to the single-particle level spacing but small 
compared to the tunneling channel band-widths. For the case of only a few modes can 
tunnel from one quantum dot to another one, the quantum point contact can be described 



by a quantum wire JTRinj. We can also use a potential barrier to adjust the transmis- 



sion coefficient. Based upon above considerations, we consider the following Hamiltonians 



15.16.1 



H = -ihv F ^2 \ dx[ip R \ T (x)d x ip Ra (x) - ip^(x)d x ip La -(x)] (1) 

H c = E c (n - |) 2 (2) 

Hi = V 2kf £^(0)^(0) + $k(0)^M0)] (3) 

where ipRcr{x) are right-moving electron operators, ipLa{ x ) are left-moving electron operators; 
E c = e 2 / (2C), eN = C g (V g R — V 9 l) is the gate voltage parameter of two quantum dots, and 
C g is a gate-to-dot capacitance, eh = (e/2) J2a Iq° d x [4 , Ra( x )i J Ru( x ) ~~ i J Lo(,~ x )' l pLa(— x )] is 
the charge difference between two quantum dots; is a backward scattering potential 
which controls the transmission coefficient T. The system described by these Hamiltonians 
(0), (0) and ([|) has been extensively studied by directly using the bosonization representation 
of the electron fields ipR(L)a{ x ) in the strategy of perturbation methods. However, the 
backward scattering term is relevant in the terminology of the renormalization group, the 
perturbation method may fail. To more effectively study this problem, we adopt another 
way developed in Ref. |14| where the backward scattering term can be exactly cancelled 



by an unitary transformation, and its effect is reflected on the correlation functions of the 
electrons and the expression of the charge operator eh. We define the following new fermion 
fields 

while the bosonic representation of these fermion fields can be written in the usual way 
1311113 as 



where D is the band width of the conduction electrons. The boson fields $1(2)0- 0*0 have the 
relation with the density operators d x ^i^,)a{ x ) — 27rp 1 ( 2 ) (T (x), and satisfy the commutation 
relations [d x $ la {x), § 1<T >{y)] = i27r8 a(T >5(x-y), [d x ^ 2 a( x ), $ 2(7 '(y)] = i2irS mr 'S(x-y), where 
the density operators are defined as pi(2)a{ x ) — V'^wO^V'i^o-t^) • The Hamiltonians (|Ip 
and (H) can be written in the bosonic representation of the fermion fields ^1(2)0- (x) as 

#0 = ^ E / ^[(«9 X .$ 1(7 ) 2 + (<9 X .$ 2(J ) 2 ] (6) 

Hi = ^T,\^Ux)-d x ^ 2a {x)]\ x=0 (7) 

where 5 = arctan(V2fc F /(^f_F)) is phase shift induced by the backward scattering potential. 
According to Eq.(|J), the charge operator eh can be written as 

1 f°° 

" = o 12 / daJ^i^)^^) + ^(^O^Mx) - if}t a (-x)ifa,r{-x) - il)tA- x )^{- x )] (8) 

O 

which only includes the cross terms of the fermion fields ipi a (x) and ip2a( x )- 

In order to cancel the 5-term in (0), we define the following unitary transformation 

f/ = exp{^E[ < fi CT (0)-$ 2 .(0)]} (9) 
Under this unitary transformation, we can obtain the following relations 

U + (H + H I )U = H 
H c = U + H C U = E c [U + nU - y] 2 (10) 

It is worth notice that the backward scattering term completely disappears, its effect is 
only reflected on the charge operator eh for this problem. It has not any influence on the 
correlation functions of the fermion fields ipi(2)a{x) because we do not consider the electron- 
electron interactions in the quantum wire. The calculation of U + nU is very simple. By 
using the following formulae 



$ 1(J ' (y)] = ivr5 ff(7 /S5(n(x - y), [$ 2(J (x), $ 2(T ' (j/)] = i^S aa/ sgn(x - y) 

6 6 
U + ® la (x)U = $i CT (x) + -sgn(x), U + $ 2a (x)U = $ 2(J (x) - -sgn(x) (11) 

and taking the following gauge transformations 

il>i„{x) = $ lc {x)e i6 \ ii2a{x)=i> 2a {xy e \ 6,-02 = 5 (12) 

we can easily obtain the following expression of the charge operator eU + hU 

U + hU = n+ CQS(2 f ^ - 1 £ / dx[^(x)^ 2ff (x) + {x)$ 1(T {x)\ 



— i- 



2 

.sin(25) 



5^ f dx[ipf a {x)^2a{x) - ip^(x)ip ll7 (x)] (13) 



It is noted that although the backward scattering term is eliminated by the unitary trans- 
formation U (0), the charge operator eU + nlI becomes complex. To simplify it, we re-define 
the right- and left-moving electronic fields 

1 - - - 1 - 

fl>JRff(x) = —^[lpla(x) + 1p2a(x)], ^La(x) = -j=[lp 1(7 (-x) - 1p 2 a(-x)) (14) 

while their bosonic representations reads ^r(l)u{x) = ( 5^7 ) ex P{~^R(L)a(x)}. These 
new right- and left-moving electronic fields are different from original ones due to we have 
taken the unitary and gauge transformations of the fermion fields ipi(2)a{x)- In terms of 
these new electronic fields, the charge operator eU + hU can be rewritten as 

U + hU = n 5 + C ° S(2 / } + 1 E(^(0) " <W0)) 



sin(25) 
— 1 — - — 



roo _ _ _ _ 

dx[lpt a (-x)lpRa(x) - ^ a (x)tp La (-x)] (15) 

„ Jo 



where n s = (l/(47r)) E ct [cos(25)4> Rr7 (00) + ^(-00) - 0l ct (oo) - cos(2<5)0 Lo .(-oo)]. If we 
take the value of the phase shift 5 as 5 C = ±vr/2, we can obtain the following relation 

U + nU\ 5 c =n 5c = \ ]T(iV iff - N Ra ) (16) 

where Nr(l)o- = are the right- and left-moving electron numbers, 

respectively. Therefore, the value of the phase shift 5 C = ±7r/2 corresponds to strong 



coupling critical point of the system induced by the backward scattering potential. At this 
strong coupling critical point, the transmission coefficient T of the electrons from the left 
quantum dot to the right one (or vice versa) is zero, i.e., the electrons are completely reflected 
on the impurity site x — 0. These two quantum dots are completely separated, therefore, 
we can observe a series of the double degenerate Coulomb blockade conductance peaks at 
this strong coupling critical point. 

We now study the low temperature behavior of the tunneling conductance and the frac- 
tional Coulomb blockade conductance peak splitting away from this strong coupling critical 
point (5 C = ±7r/2). Because the directly hopping of the electrons from the right of the 
impurity to its left or vice versa is very weak, we can approximately replace the quantity 



oc 



dx[^i a (-x)il) RrT {x) - i)l a (x)i) La (-x)} 
by 

a[^(0)^(0) - ^(0)^{ 



where a is a small constant. For simplicity, we define the following new boson fields 

In terms of these new boson fields, the charge operator eU + nU is 

cos(25) + 1 . , aDsm(25) ... /T ... . , 

U+nU = n 5 + V + — sin(<l>_ c (0)) cos($_ s (0)) (17) 

and the Hamiltonian (||) reads 

H = ^ J dx[{d x ^ +c {x)f + (d x $ +s (x)) 2 + (d x ^ c (x)) 2 + (d x ^s(x)) 2 } (18) 

However, due to the quantity a sin(2<5) is very small, the last term in flTT| ) can be perturba- 
tively treated, the electrostatic Hamiltonian (0) can be written as by using Eq. fll7D 



He = £c[ COS(2 / ) + V c (0) - AAf + Ec ( aD f m r sin 2 ($_ c (0)) cos 2 ($_ s (0)) (19) 

where AiV = — + N/2. It is worth noting that for the case of weak backward scattering 
5 ~ 0, i.e., the transmission coefficient T ~ 1, the low temperature behavior of the system 
is different from that for the strong backward scattering. We first consider the case of the 
weak backward scattering 5 ~ 0. Because the factor (cos(2<5) + l)/(27r) is finite, the boson 
field $_ c (0) has a non-zero expectation value 

< *_„«» >= AN (20) 

which suppresses its low energy excitations. Therefore, we can safely integrate out the boson 
field <3>_ c (0) and obtain the following effective Hamiltonian 

-/ / aZ^sin(25) N9r Ane^Ec^ , , , . . ... , , 

Hc = Ec{ 2nhv F ) [1 ~ (_ D~^ ] C ° S(2 < $ - c(0) >)][1 + cos ( 2$ -( ))] ( 21 ) 

where E c = E c (l+cos(25)) 2 / (2tt) 2 , and 7 ~ 0.577 is Euler constant. However, in the fermion- 
ization representation, cos(2$_ s (0)) can be written as: cos(2$_ s (0)) ~ (Q + )i[> + (Q~) + 
^(0-)^(0 + ), where ± = ±p, p is an infinitesimal quantity, which produces a gap in the 
spectrum of fermion ip(0). Therefore, the cos(2$_ s (0))-term in fl2"ip is irrelevant in the 
terminology of the renormalization group. 

For the case of the strong backward scattering 5 ~ ±tt/2, the electrostatic energy signif- 
icantly depends upon the quantity AiV. If AiV is zero, AiV = 0, after integrating out the 
boson field $_ c (0), we can obtain the following electrostatic Hamiltonian 
- T / _ ,aD sin(25) , 9r „ Ai\e 1 E c . 9l r „ 

For the case of AiV 7^ 0, because of |Max($_ c (0))| = 7T, if the phase shift 5 satisfies the 
relation: \S\ < 5 , where 5 is defined as: 1 + cos(25 ) = 2 AN, after integrating out the 
boson field $_ c (0), we can obtain the same Hamiltonian as that in (|21|). If the phase shift 
satisfies the relation: 5 < 1^1 < 7r /2, we can obtain the following effective Hamiltonian 
s , cos^J + l <<S ,_ M> _ AN? 

^ ,aD sin(2<5) . 9r „ . 47re 7 ^ c . 9l r „ 
+ 2?rfi J T [l - (— ^) 2 ][1 + cos(2$_ s (0))] (23) 
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where < $_ c (0) >= ti for AN > 0, and < $_ c (0) >= -7r for AiV < 0. It is worth 
notice that the present method used to get the final effective Hamiltonians in (pTD, (F22]) and 



(p3|) is completely different from that in Refs. (15| and ||21[| . After performing the unitary 
transformation (9), all information of the conduction electron scattering on the impurity 
is incorporated into the electrostatic Hamiltonian (19) which only contains higher order 
terms of the boson fields <3>_ c ( s )(0). However, in Refs. |15| and [|2lj, the authors directly 



take the approximation on the backward scattering potential which leads to a lower order 
terms of the boson fields $_ c ( s )(0). This difference may make us obtain a different transport 
behavior from that in Ref. |2l[] in the low temperature limit. On the other hands, in terms of 



the terminology of the renormalization group, the former only gives some irrelevant terms, 
but the latter gives some relevant terms. Based upon above discussions, for the case of 
< \5\ < 5 , we can obtain the following ground state energy of the system 

Tivf 

For the case of 5q < \5\ < tt/2, the ground state energy of the system reads 

E = E + E C [\AN\ - i(cos(25) + l)] 2 

- ge( ggg§iMgj) )a (25) 

Now we can obtain the fractional Coulomb blockade conductance peak splitting jnj. For 
the case of < \8\ < 5 , it can be written as 

/ = [E(AN = i) - E{AN = 0)\/E(AN = \,S' = ±|) 

= 2ae^sin( M ) l + o^ o 

Tiv F 2it 1 + cos(2d) 

In the case of <5o < |5| < tt/2, it can be written as 

/ = 4[i-i(l + cos(25))] 2 (27) 

It is noted that the fractional Coulomb blockade conductance peak splitting significantly 
depends upon the phase shift 5, i.e., the transmission coeficient T. However, in the case of 
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the weak backward scattering < \8\ < 5q, it also denpends upon the parameter a, the 
Fermi velocity Vp and the electrostatic Coulomb energy E c = e 2 /(2C). In the case of the 
strong backward scattering 5q < \5\ < rr/2, the fractional Coulomb blockade conductance 
peak splitting / only depends upon the phase shift 5, therefore, it is universal. 

We now study the low temperature behavior of the tunneling conductance. To this end, 
we define the following current operator 

I = ed t (U + nU) = e 1 + C o ° s(2 ^ ^ c (0) 
_ aeDwi(2S) C os($_ c (0)) cos($_ s (0)) - d t ^ s (0) sin($_ c (0)) sin($_ s (0))] (28) 

TTllVp 

Using the Kubo formula of the conductance, we can easily obtain the following tunneling 
conductance in the low energy and low temperature limit 

G(T) = G (5)[l + AT( S -^-) 2 cos(2 < $_ c (0) >)] (29) 

7T 

where Go(5) = e 2 (l + cos(25)) 2 / (4^h), T is the temperature, and A is a constant. It is noted 
that the tunneling conductance has a linear temperature dependence, and at the strong 
coupling critical point 5 C = ±7r/2, it is equal to zero. It is necessary to mention that as 
the frequency to and the temperature T tend to zero, the system approaches to the strong 
coupling critical point delta = ±7r/2 because the backward scattering term is relevant, the 
renormalized backward scattering potential V2k F tends to infinity in the low energy limit 
22|. 

We define the following dimensionless tunneling conductance 

g = ^ = \(l + cos(26)Y (30) 

where Go — e 2 / (2ttH) is a unit quantum conductance of each tunneling channel. From 
Equs.(ppT) and (0), we can obtain the following relation between the dimensionless tunneling 
conductance g and the fractional Coulomb blockade peak splitting / 

( 8V2BgV 2 (l - M(l - cos(^=)), < \S\ < S 

/= 1 V2 ^ (31) 

( 4(- - f-) 2 , 5 < \S\ < vr/2 



where B = [ae 1 E c / '{t^Tlvf)} 2 '■ It qualitatively agrees with the experimental data in Ref. 

Because the factor a depends upon the structure of the quantum point contact, 
and the electrostatic energy E c is determined by the structure of the quantum dot, the 
fractional Coulomb blockade peak splitting is not universal in the weak backward scattering 
< \S\ < Sq. However, in the strong backward scattering 5q < \S\ < tt/2, it is universal. 

In summary, by using the bosonization method and the unitary transformation, we have 
shown that the tunneling conductance between two quantum dots has the linear tempera- 
ture dependence in the low energy and low temperature limit, and given a general relation 
between the dimensionless tunneling conductance g and the fractional Coulomb blockade 
peak splitting / which is valid both for the weak and strong backward scattering of the 
electrons on the impurity. Our treatment of this kind of system is simple and effective, 
and can be used to study other similar problems because it can exactly treat the backward 
scattering potential which is a relevant term in the low energy region, the usual perturbation 
expansion methods of this term may fail because the high orders are divergent in the low 
energy limit. 

We are very grateful to Prof. Peter Fulde for his encouragement. 
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